تحلیل روند تغییرفاز ناشی ازتغییر ساختار کریستالی آلیاژ های حافظه دارشکلی تحت بارهای مکانیکی یا حرارتی بسیار پیچیده است. با توجه به درگیری ترمومکانیکی موجود در آلیاژ های حافظه دار شکلی، در صورت سریع بودن نرخ بارگذاری، تولید/جذب حرارت طی تغییرفاز رفت/برگشت منجر به تغییرات دمایی ماده می شود که بر رفتار مکانیکی آن تاثیرگذار است. در این مقاله الگوریتمی مبتنی بر روش اجزاء محدود جهت درنظر گرفتن رفتارهای پیچیده مکانیکی، حرارتی و درگیر آلیاژ های حافظه دارشکلی ارائه شده است که هر دو رفتار منحصر به فرد این مواد یعنی سوپرالاستسیته و حافظه دار شکلی را در بر می گیرد. به منظور نشان دادن کارایی و دقت الگوریتم ارائه شده، چند مساله کلیدی شبیه سازی و مورد بحث قرار گرفته است.